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O ' We present 850 /zm polarimetry from the James Clerk Maxwell Telescope to- 

ward several dense cores within the dark cloud Barnard 1 in Perseus. Significant 
polarized emission is detected from across the mapped area and is not confined 
Oh" to the locations of bright cores. This indicates the presence of aligned grains and 

hence a component of the magnetic field in the plane of the sky. Polarization 
vectors detected away from bright cores are strongly aligned at a position angle of 
~ 90° (east of north), while vectors associated with bright cores show alignments 
of varying orientations. There is no direct correlation between the polarization 

^ ■ angles measured in earlier optical polarimetry toward Perseus and the polarized 

submillimeter thermal emission. Depolarization toward high intensities is exhib- 
ited, but toward the brightest core reaches a threshold beyond which no further 
decrease in polarization percentage is measured. The polarized emission data 
from the interior envelope are compared with previously published OH Zeeman 
data to estimate the total field strength and orientation under the assumption 
of a uniform and non-uniform field component in the region. These results are 
rough estimates only due to the single independent detection of Zeeman splitting 
toward Barnard 1. The uniform field component is thus calculated to be B = 31 
fiG [±(0.52iV — 0.01E) — 0.86i] in the case where we have assumed the ratio of 
the dispersion of the line-of-sight field to the field strength to be 0.2. 



Subject headings: ISM: clouds, magnetic fields, molecules — polarization — ISM: 
individual (Barnard 1) — stars: formation — submillimeter 
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1. Introduction 

Over the past decade, observational evidence for the presence of magnetic fields in molec- 
ular clouds and their role in star formation has grown dramatically. However, detection of 
magnetic fields is not synonymous with measuring their geometry or the total field strength. 
Recent simulations have shown that the low density regimes of molecular clouds may not 
be in magnetic equilibrium (Padoan & Nordlund 1999), although in high-density cores, the 
evidence for magnetic and virial equilibrium is stronger (Myers & Goodman 1988; Crutcher 
1999; Basu 2000). While models can help interpret data, it is still very rare for evidence of 
Zeeman line splitting (which traces the line-of-sight field, Bi os ) and thermal dust polarization 
(which traces the plane-of-sky field, B pos ) both to exist toward regions of similar density in 
a single cloud. 

The Perseus molecular cloud complex is one of the closest star-forming regions to the 
Sun. Its distance is the subject of some debate, but it is thought to be associated with the 
Per OB2 association at a distance of 334 ± 12 pc (Borgman & Blaauw 1964). However, the 
complex is likely in front of the OB association (Lynds 1969), and Cernicharo et al. (1985) 
suggest that there are in fact two clouds along the line of sight toward the complex, the 
second at a distance more comparable to Taurus at 200 pc. In this paper, we adopt a 
distance of 330 pc for the Perseus association. CO emission reveals that the complex is 
elongated, extending over 55 pc along its major axis (at ~ 60° east of north), but just 15 
pc along the minor axis (Sargent 1979). Along its length, six denser star-forming clouds 
(L1448, L1455, NGC 1333, Barnard 1, IC 348 and Barnard 5) are connected by low density 
molecular gas of n ~ 10 2 cm -3 (Bachiller & Cernicharo 1986). 

The Barnard 1 (Bl) cloud has been observed in many molecular transitions and modeled 
as a multi-phase cloud with a thin outer envelope, denser inner envelope and a central dense 
core. Bachiller & Cernicharo (1984) observed the cloud in several isotopes of CO, HCO + 
and a single transition from NH 3 and determined that temperatures are higher toward the 
outer edges of the cloud, indicating primarily external heating. Three optically visible young 
stellar objects, LkHa 327, LkHa 328, and LZK 21 are associated with the cloud, two of which 
show IRAS emission. Three additional IRAS sources are undetected optically; the presence 
of these sources indicates some recent star formation in Bl. In the center of the cloud, 
Bachiller et al. (1990) observe CS J = 1 — emission from dense gas over a region ~ 2 pc x 
5 pc with an accompanying mass of 1200 M . NH 3 (1,1) and (2,2) emission toward the CS 
"main core" reveal substantial substructure within the gas, showing evidence for two or three 
condensations (Bachiller et al. 1990). IRAS 03301+3057 lies at the center of the main core, 
but does not coincide directly with any of the ammonia peaks; it is located approximately 1' 
north of the south-western ammonia peak (the "southern clump" of Bachiller et al. (1990)). 
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High resolution observations of H 13 CO + by Hirano et al. (1999) reveal a strong peak 
(denoted Bl-b) at the position of the south-east ammonia emission detected by Bachiller et 
al. (1990). As part of the same study, continuum emission at 850 fim 350 fim (from the 
JCMT and CSO, respectively) and 3 mm (from the Nobeyama Millimeter Array) clearly 
identify two high density dust cores within the single molecular clump. Based on spectral 
energy distributions, Hirano et al. (1999) conclude that these objects are both extremely 
young protostars in the Class phase (Andre et al. 1993). The masses of the central 
objects, called Bl-bN and Bl-bS, are estimated to be no greater than 7x 10~ 2 M , indicating 
extremely young ages of less than 2 x 10 4 yr for both sources. No outflows have been detected 
from either source. 

One powerful outflow, associated with IRAS 03301+3057 (Bl-IRS), has been identified 
by Bachiller et al. (1990) in CO (J = 2-1). It is confined to a region of 40" (about 0.07 pc). 
The dynamical time estimated from the outflow is 10 3 to 10 4 years. Hirano et al. (1997) 
measured the small scale structure of this CO outflow and estimate that the driving source 
is very young and is observed in a pole-on configuration. 

The masses of the YSOs seen optically and with IRAS range from 0.2 to 3 M Q (Bachiller 
et al. 1990). Given these low stellar masses, the stellar to gas mass ratio in Bl is ~ 0.5%, 
negligibly small even compared with Taurus where the star formation efficiency has been 
recently estimated at 6% (Onishi et al. 1998). The observed rotation velocities within Bl 
are insufficient to support the cloud against collapse by a factor of ~ 8 (Bachiller et al. 1990). 
The ages of embedded but optically visible objects LkHa 327 and LkHa 328 are between 
4-6 xlO 6 yr (Cohen & Kuhi 1979). Based on this, Bachiller et al. (1990) conclude that 
a mechanism must be providing substantial support to the Bl cloud. The mechanism is 
generally attributed to a magnetic field. 

Polarization of background starlight from the Perseus cloud due to selective absorption 
from dust grains within the complex was measured by Goodman et al. (1990a), who find 
that the distribution of polarization position angles is bimodal, with weaker vectors (less po- 
larized) aligned along the cloud's projected major axis and stronger vectors (more polarized) 
lying roughly perpendicular to the first population. Goodman et al. (1990a) hypothesized 
that two clouds of differing magnetic field orientations could be superimposed along the line 
of sight. A prior argument for a second gas cloud along the line of sight to Perseus and Bl 
at a distance of 200 pc was presented by Cernicharo et al. (1985). 

The Bl cloud has been surveyed for evidence of Zeeman splitting in dense OH gas more 
extensively than any other dark cloud. Lang & Willson (1979) estimated a 3a limit of 90 
/xG toward LkHa 327, located approximately 4' away from Bl's strong molecular peak. Bl 
was chosen by Goodman et al. (1989) as a strong candidate for magnetic field detections 
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due to its atypically high non-thermal linewidth components, and a field strength of —27 ±4 
/j,G was measured toward the position of the bright molecular core coincident with IRAS 
03301+3057. (The negative sign indicates that the field is oriented toward the observer.) A 
survey of 12 dark clouds for evidence of Zeeman splitting yielded only one solid detection - 
toward Bl - with the 140 ft. Green Bank Telescope (Crutcher et al. 1993). In an observation 
toward the source IRAS 03301+3057, a field strength of —19 ± 4 jiG was measured, which 
is consistent with the Arecibo value when beam dilution is taken into account. 

In order to supplement the Zeeman data toward the dense molecular gas of Bl, we have 
measured polarized emission at 850 /zm from dust toward the "main core" of Bl as identified 
in CS and NH 3 , both tracers of high column densities. Emission from aligned, spinning 
dust grains is anisotropic and hence polarized. Unfortunately, polarization data reveal no 
direct information about the field strength, since the degree of polarization is dependent 
on other factors such as grain shape, composition and degree of alignment. The degree 
of polarization is in essence a measure of how effectively the grains have been "sped up" 
(Hildebrand et al. 2000). However, even though the grain spin is induced by mechanisms 
other than the magnetic field, such as the radiation field (Draine & Weingartner 1996) or 
the production of H 2 on the grain surface (Purcell 1979), the magnetic field is expected to 
provide the alignment. Because of this, continuum polarization data are the principal means 
of probing the geometry of the magnetic field. The very sensitive Submillimetre Common- 
User Bolometer Array (SCUBA) detector now permits the observation of polarized emission 
from the ambient cloud surrounding dense cores. 

This paper is the fifth in a series to examine the magnetic field geometries in star- 
forming molecular clouds using polarized emission at 850 /im. Barnard 1 is the first dark 
cloud we have observed, and these data are the first emission polarimetry toward this region. 
The observations and data reduction techniques are described in § 2. The polarization data 
are analyzed in § 3. We discuss the possible interpretations of these data and calculate 
an estimated three dimensional structure for a uniform field component of Bl in § 4. Our 
findings are summarized in § 5. 

2. Observations and Data Reduction 

We have used the UK/Japan polarimeter with the SCUBA detector at the James Clerk 
Maxwell Telescope 1 , to map polarized thermal emission from dust at 850 fim toward a dense 



1 The JCMT is operated by the Joint Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the UK, the Netherlands Organization for Scientific Research, and the National Research 
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region of the Bl dark cloud. The observations were taken from 11 to 13 October 1999. 
The polarizer and general reduction techniques are described in Greaves et al. (2000) and 
Greaves et al. (2002). To generate a polarization map, a 16-point jiggle map was made at 
each of 16 different half-waveplate positions. After each 12s integration, the half-waveplate 
was rotated through 22.5° and the mapping repeated. The data were flat-fielded, corrected 
for extinction and dual-beam corrected using the standard SCUBA software. Estimation of 
systematic errors due to chopping and sky subtraction can be found in (Matthews, Wilson 
& Fiege (2001), hereafter Paper II). Unfortunately, there are no large scale 850 /iin maps 
of the Bl cloud in the literature. This made the identification of chop angles and sky 
removal candidate bolometers more difficult. The chop angles and throws for each field 
center observed are summarized in Table 1. The level of precipitable water vapor was very 
stable over the course of the observations. The estimates of r(225 GHz) from the CSO 
ranged from 0.055 to 0.075 over the observations, with 96% in the range 0.060 to 0.070. 

These data were corrected for an error in the SCUBA clock which placed incorrect LST 
times in the data headers during the period from July 1999 to May 2000. This error did not 
affect the telescope's acquisition or tracking, but affects data reduction since the elevation 
and sky rotation are calculated from the LST times in the data headers. The magnitude of 
this error over time can be evaluated and then corrected retroactively as described on the 
JCMT website. The error in timing after this adjustment is ±10s. The data were reduced 
using the Starlink software packages POLPACK and CURSA, designed specifically to include 
polarization data obtained with bolometric arrays. 

After extinction correction, noisy bolometers identified for each night's data were flagged 
and removed from the data. Between 3 and 5 bolometers were removed per night. Prior to 
sky subtraction, images were made to examine the flux in each bolometer, since bolometers 
used for sky subtraction should not have negative values (produced if one has chopped onto a 
location with significant flux, for example). The data were sky subtracted using bolometers 
with mean values close to zero. Between 1 and 4 bolometers were used to subtract the 
sky. The methods of sky subtraction are discussed in detail in Jenness, Lightfoot & Holland 
(1998). Finally, the instrumental polarizations (IPs) were removed from each bolometer. All 
the data sets were then combined to produce maps of three Stokes' parameters (/, Q, and 
U), which were then combined to yield the polarization percentage and polarization position 
angle according to the following relations: 

V = ^ Q2 j +U2 ; 6 = X - arctan([//Q). 



Council of Canada. 
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The uncertainties in each of these quantities are given by: 

dp = p~ l ^[dQ 2 Q 2 + dU 2 U 2 }; d6 = 28.6° ja v 

where a p is the signal-to-noise in p, or p/dp. 

A bias exists which tends to increase the p value, even when Q and U are consistent 
with p — 0, because p is forced to be positive. The polarization percentages were debiased 
according to the expression: 

Pdb = vV - dp 2 . 

Future references to polarization percentage, or p, refer to the debiased value. 

Absolute calibration is not part of the standard reduction of polarization data since 
the percentage polarization is a relative quantity. However, from our Stokes' / map, we 
can estimate fluxes by using a reasonable flux conversion factor for 850 /im SCUBA data. 
This quantity is dependent on the chop throw used, and for a throw of 120", the standard 
flux conversion factor is 219 ± 21 Jy beam^V" 1 according to the JCMT website. However, 
the flux conversion factor is at least a factor of two greater for polarization data due to 
the presence of the analyzer which retards half the incoming flux on average. In practice, 
the flux conversion factor is ~ 2.2 times the standard value due to imperfect transmission 
through the waveplate (J. Greaves, 2002, private communication). Hence, we have scaled our 
Stokes' / map by 480 Jy beam _1 V _1 . The associated uncertainty in this factor introduces 
an uncertainty of 10% into the resultant fluxes. 

Before filtering the data to select reliable polarization vectors, it was necessary to esti- 
mate the effects of sidelobe polarization in the position of the main beam. This is a measure 
of the minimal believable polarization, p cr u, given the potential for sources in sidelobes 
to produce artificial polarization signals in the central region of the map (see Greaves et 
al. (2002)). For our worst case scenario in Bl, the flux contributed at approximately 68" 
from the map center is 16 times that at the center. Examination of polarization maps of 
Saturn (which has only a small intrinsic polarization ~ 0.6%, with a minimum of ~ 0.2%) 
of 13 October 1999 reveal that the relative mean power 68" from the main beam center is 
0.0067 compared to the main beam itself. The mean polarization percentage at this position 
on the SCUBA field is 4.3%, which is a measure of the instrumental polarization. The p crit 
value is given by: 

Peru > 2 x 4.3% x 0.0067 x 16 

(see Greaves et al. (2002)) which gives a minimum threshold polarization of 0.92%. Taking 
into account the intrinsic polarization of Saturn, this leaves ~ 0.3—0.7% arising from sidelobe 
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polarization. We have thus selected vectors for which polarization percentage, p > 1% as 
reliable data. The vectors selected also have an uncertainty in polarization percentage, 
dp < 1%, and signal-to-noise in polarization percentage a p > 3. To minimize the systematic 
effects arising from the possibility that we have chopped onto a region of polarized emission, 
vectors are selected only if they are coincident with Stokes' / > 20% of the faintest peak in 
our map. As discussed in Paper II, if the reference position has a flux level ~ 10% that of 
the source peak and polarized to the same degree, in the final map the position angle at 20% 
the peak would be offset from the correct value by < 10°, while the p value is incorrect by 
at most a factor of two. For the brighter peaks, the effects would be considerably reduced. 



3. 850 |ira Polarization Data 

Figure 1 illustrates the polarization pattern detected across the Bl "main core" region 
as identified in CS and NH 3 by Bachiller et al. (1990). The polarization data are plotted on 
a colored greyscale Stokes' I map estimated by summing together the fluxes detected at all 
waveplate positions. Table 2 contains the data (with p > 1%) in tabular form. Four peaks 
are distinguishable, these are labelled Bl-a to Bl-d. Bl-a and Bl-b (N and S) follow the 
classification of Hirano et al. (1999) as identified in H 13 CO + and 850 /im SCUBA emission. 
The presence of two sources within Bl-b was confirmed by 3.0 mm observations with the 
Nobeyama interferometer (Hirano et al. 1999). NH3 (1,1) and (2,2) emission was observed 
from peaks corresponding to Bl-a, Bl-b and Bl-c by Bachiller et al. (1990). The Bl-d 
850 /im peak lies approximately 1' south of the Bl-a molecular peak (Bachiller et al. 1990) 
which is likely associated with IRAS 03301+3057 (marked by a blue cross on Figure 1). No 
ammonia emission is concentrated at the Bl-d position (Bachiller et al. 1990), although 
a very low signal-to-noise peak exists near this position in the H 13 CO + map of Hirano et 
al. (1999). 



3.1. Polarization Position Angles 

Polarized emission is detected both on the bright cores and in regions of lower column 
density between them. The degree of alignment across the region is evidence for the presence 
of ordered magnetic fields within the main core of the Bl cloud. The data in high intensity 
regions have been binned to 6" sampling, while data in fainter regions are binned to 12" 
to improve the signal-to-noise ratio. The distributions of polarizations associated with faint 
emission and bright emission are plotted separately on Figure 2. 
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The distribution for faint / (dashed line) is approximately Gaussian. A fit to these data 
yields a mean of 91.3° with a distribution width of 19.0°. Polarization position angles are 
measured such that values increase east of north. A goodness of fit measure to the data yields 
xfed = 0-6- The statistical mean of the distribution is 88.3° (east of north) with a standard 
deviation of 27.7°. The distribution of vectors in regions of bright emission, however, cannot 
be fit effectively by a Gaussian (or even a series of Gaussians). The solid line of Figure 2 
shows several peaks, each of which corresponds roughly to one of the bright peaks. We have 
indicated the peak sampled on the distribution. 

Therefore, the polarization pattern in the ambient cloud material around the cores is 
defined by a mean polarization direction where the vectors are distributed about 90° (east 
of north), while the cores each exhibit different mean position angles. The core Bl-b shows 
systematic variation in position angle. The northern core has < 9 >~ 65°, while the southern 
core exhibits < 9 >~ 120°. The Bl-c core has < 9 >~ 35°, and the Bl-d core peaks around 
90° east of north, in alignment with the fainter material in the cloud. This could indicate 
that this core is not strongly differentiated from the ambient cloud yet. 

On Figure 2, we have also plotted the position angles associated with the two peaks 
observed in the optical absorption polarization data of Goodman et al. (1990a), taking into 
account the 90° offset expected for emission polarimetry. These positions (55° and 161°) 
are not coincident with any of the peaks in the 850 /zm distribution, either at high or low 
intensity. 

3.1.1. Correlations between Adjacent Vectors 

To better examine the changes in the nature of the polarization data across the map, 
we have compared each vector to its closest eight neighbours, calculating the differences in 
polarization percentages and position angles for each pair of values. All the data of Figure 
1 were used, except those values for which p < 1% (shown in red). Next, the results were 
smoothed onto 12" grid, by calculating the mean changes in polarization percentage and 
position angle. The resulting map is shown in vector form in Figure 3, where the vector 
magnitude is the mean change in polarization percentage in a grid unit and the vector angle 
is the mean change in orientation. No significant change in orientation is indicated by a 
vector at 0° east of north. 

Total intensity contours illustrate the positions of the cores on Figure 3. Only one peak 
is distinguishable in Bl-b, but it is clear this core is elongated. As expected, the changes 
within the cores are relatively small. The position angles (and even polarization percentages) 
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are consistent with relatively little change. Based on the histogram of the polarization data 
toward fainter regions shown in Figure 2, the vectors in this regime were also expected to be 
well aligned, and the data of Figure 3 show that this is indeed the case. The large variations 
in adjacent grid units are confined to the edge of the map and the boundaries between the 
Bl-b and Bl-c cores and the fainter material. 



3.2. Continuum Fluxes Toward Cores in Bl 

Of the four dense cores detected in our polarization map, two have not been observed 
previously in continuum emission. These are Bl-c and Bl-d, although the latter may have 
been confused with Bl-a in large beams (i.e., IRAS) particularly if they are at similar evolu- 
tionary stages. The brightest source in our map is Bl-c, as revealed by the Stokes' I contours 
of Figure 3. We do not detect both peaks toward Bl-b although Hirano et al. (1999) do in 
an earlier SCUBA map; their interferometic observation with the Nobeyama Millimeter Ar- 
ray clearly resolves two peaks toward this source. The slightly enhanced emission coincident 
with IRAS 03301+3057 is the faintest distinguishable peak at 850 /im. Table 3 summarizes 
the peak fluxes toward each of the four cores (for Bl-b, the peak flux is that of Bl-bS) and 
positions of these peaks. 

For a cloud at a distance of 330 pc and assuming /tsso/xm = 0.01 cm 2 g" 1 , a mean 
molecular weight of 2.33 and a dust temperature of 20 K, we find that the column density 
N(H2) = 10 23 S < 85o[Jy beam -1 ] cm -2 . The Bl-c peak has a flux of 3 Jy beam -1 , which implies 
N(H 2 ) = 3 x 10 23 cm -2 , or 300 magnitudes of visual extinction. Toward Bl-bS, we estimate 
N(H 2 ) = 2.5 x 10 23 cm -2 , which is within 50% of the estimate of Hirano et al. (1999). 
Assuming a core depth comparable to the FWHM of 30" for Bl-c yields a volume density of 
2 x 10 6 cm" 3 which is typical of prestellar and protostellar core densities. 

3.3. Depolarization in Barnard 1 

Figure 3 suggests that changes in polarization percentage are small within the cores of 
Bl. The statistical means of the low intensity and high intensity vector populations are 4.5% 
and 2.6% in 81 and 74 values respectively. The standard deviations in these populations are 
2.3% and 1.4%. In this case, the depolarization effect, declining polarization percentage 
with increasing intensity, may be weak within parts of Bl. The easiest way of examining the 
depolarization effect is to plot p versus / for all vectors on the polarization map. 

In Figure 4, we plot the data for the Bl region as presented in Figure 1 excluding only 
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those data values with p < 1% (plotted in red). The data exist in two populations, where the 
data at low intensities are binned to 12" (shown as crosses) and the data at high intensities 
are binned to 6" (shown as circles). Although these plots are shown in log-log space, the fits 
to the data were done to profiles of p versus / by minimizing x 2 - This is a more effective 
treatment of the uncertainties since those for low values of p are exaggerated in log space. 
The fits to these two populations produce completely consistent slopes, indicating that both 
can be characterized by power laws of the form: p = AP with an index of 7 « —0.8. At 
high values of / however, there is a slight thresholding of polarization percentage. 

To better illustrate this, in Figure 5 we show the same style plot for three cores: Bl-b 
(north and south combined), Bl-c and Bl-d. Vectors toward the cores Bl-b and Bl-c exhibit 
higher values of p than expected given the declining trend below 1 Jy beam -1 . In fact, the 
distribution of p versus / flattens at high intensities. In the case of Bl-b, this flattening 
could be the result of our lack of sensitivity to values of p < 1%. However, for the Bl-c 
core, this constraint removes only a single vector (which has a value of 0.99 ±0.25%). Thus, 
in Bl-c, the depolarization effect does not follow the usual trend of declining polarization 
percentage as intensity increases. (The Bl-d core does exhibit depolarization to its peak, 
which is significantly lower in intensity than Bl-b and Bl-c.) The Bl-c core is well sampled 
(with only one vector missing) and definitely exhibits a flat dependence of p on / down to 
30% of that core's peak. The omitted vector (shown in red on Figure 1) does not correspond 
to the intensity peak of the core, but is associated with an intensity just 2/3 of the peak 
observed. To our knowledge, this is the first case of a bright core which does not exhibit 
depolarization over its whole observed intensity range. The observed threshold is not the 
effect of optical depth. We estimate that the optical depth at the Bl-c peak is ~ 0.015 which 
is << I. 



4. Discussion 

4.1. Interpreting the Polarization Pattern 

Optical polarimetry using absorption of light from background stars was used to probe 
the magnetic field structure through dust at low extinctions in the Perseus cloud complex 
by Goodman et al. (1990a). They found a bimodal distribution of polarization vector 
orientations toward the complex such that the vectors lie roughly parallel and perpendicular 
to the major axis of ~ 60°. They fit two Gaussians to their data set with means 71° 
and 145° east of north and 1 a dispersions of 12° and 8°, respectively. Since there is no 
spatial distinction between the two populations and evidence in observations of molecular 
gas that two different clouds could lie along the same line of sight, they conclude that the 
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two polarization populations are representative of two distinct clouds at different distances. 
The foreground cloud was predicted to have a low extinction (A v < 1 mag). Bimodal 
distributions in polarization vectors had been previously noted by studies toward the Perseus 
clouds NGC 1333 (Turnshek et al. (1980); Vrba, Strom & Strom (1976)) and Barnard 5 
(Joshi et al. 1985). 

In the process of generating an evolutionary model for the Bl cloud, Crutcher et 
al. (1994) adopt a mean plane-of-sky field direction along the minor axis of the Perseus 
(and hence Bl) cloud which corresponds to the polarization distribution centered on 145° 
east of north. In this case, the Goodman et al. (1990a) vectors centered on 71° east of north 
would most likely be associated with the foreground cloud, at 200 pc distance (Cernicharo 
et al. 1985). If the interior of the dense Bl main core is threaded with the same field geom- 
etry as measured on the periphery of the Bl/Perseus cloud, then a field of mean direction 
145 ±12° should produce polarized emission from dust at a position angle of 55 ±12°. Figure 
2 demonstrates that there is no peak in the polarization data at 850 /zm at 55°, although 
polarization vectors of Bl-c and Bl-bN fall roughly in the range of the optical polarization 
dispersion value. The second distribution measured by Goodman et al. (1990a) peaked at 
71°, which corresponds to an emission polarization angle of 161°. Thus, there is no compo- 
nent of position angle in our data set which corresponds directly to either population of the 
optical polarization data. Hence, based on our emission polarimetry, which is believed to 
arise only in the regions of dense, cold dust, we cannot conclude which optical polarization 
direction is more likely associated with the Perseus complex. 



4-1.1. Complete Depolarization of the Cores? 

The polarization vectors across the Bl-d peak are aligned with the faint emission polar- 
ization angles of 90°, but the brighter peaks, Bl-b and Bl-c, both exhibit different position 
angles. Given the previous observations of depolarization in bright peaks and the polariza- 
tion plateau across these two cores, it would be tempting to think that the cores themselves 
are completely depolarized and that one might thus measure the polarization toward a dif- 
ferent cloud along the line of sight at their positions. There are several reasons why this is 
unlikely to be the case. 

In most well-sampled cores, depolarization is a non-linear function of intensity, so the 
polarization percentage rarely reaches zero (e.g. Henning et al. (2001); Ward-Thompson et 
al. (2000); Greaves et al. (1999); Greaves, Murray & Holland (1994); Minchin & Murray 
(1994)), and if it does become unmeasurably small, as it does at places in Bl-b, then it does 
so at the highest intensities. Therefore, one might expect to see a varying position angle 
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across cores that reflects the varying contribution of the core's polarized intensity to the 
vector sum with the second cloud, revealing the second cloud's position angle only toward 
the peak of the core. Furthermore, even in the case where the depolarization might atypically 
be complete across a core, this would effectively create a steeper than usual depolarization 
effect rather than a flat distribution of p with / as we particularly observe toward Bl-c (see 
Figure 5). The reason is that if the polarized emission arose in a cloud other than Bl and 
were unassociated with the Bl-c core (which we know to be in Bl through its associated 
molecular emission), the increasing contribution to the total intensity from the core (despite 
its contribution of zero polarized intensity) should produce a variation in p with I which 
varies exactly as the increase in intensity across the core. The only way this could be 
avoided would be if, as the intensity across the Bl-c core increased, so did the polarized 
intensity in the second cloud. These increases would have to match each other exactly to 
produce the flat pattern we see. We dismiss this scenario as far-fetched. 

In Bl, we see not just a single core which exhibits a different position angle orientation 
than the ambient cloud, but instead there are two such cores, where Bl-b includes systematic 
variation in 9 from north to south. In order to see the optical polarizations from a foreground 
and background cloud in absorption, Goodman et al. (1990a) point out that the extinction 
of the nearer cloud must be low (Ay < 1 mag) or have significant fluctuations in order to 
see through to the background cloud. If the extinction in the cloud is on the order of 1 mag, 
then the cloud is not self-gravitating and hence is an unlikely source of 850 /iin emission 
(and certainly 850 fxm polarized emission, which is at most on the order of 10% of the total 
flux). This leaves the possibility of a denser cloud with significant fluctuations in column 
density across its projected surface. 

The angular distance between the Bl-b and Bl-c cores is approximately 2 arcminutes. 
At a distance of 200 pc, which is the most probable location of a foreground cloud given 
previous molecular line observations, this corresponds to 0.12 pc, or about the scale of a 
self-gravitating core in a molecular cloud. Since we observe different position angles in the 
two bright cores (and a change in position angle within the Bl-b core itself), this would 
suggest that fields of varying orientations are exhibited in the foreground fluctuations in the 
nearby cloud. This implies that the fluctuations in the foreground cloud are comparable to 
or smaller than the scale of a starless core. Since the two bright cores are associated with 
NH 3 emission (and hence known to be within the Bl cloud), the fluctuations would have to 
vary on an angular scale similar to the separation of cores in Bl, despite being at a distance 
nearly a factor of two smaller. We conclude that attributing the polarization seen against 
the bright cores Bl-b and Bl-c to a foreground cloud suggests an unlikely configuration for 
a nearby cloud not detected in any dense tracers such as CS and NH 3 . 
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However, the flattening of the p versus / relation in Bl-c does suggest that this core is 
unusual in some way. Depolarization within cores is usually attributed to either changing 
grain or alignment physics with increasing density or varying (ordered or disordered) mag- 
netic field geometry. It is possible that no stellar condensation has formed at the center 
of Bl-c, in which case a field with a straight geometry may thread this core. Also, if no 
outflow is present, then theory suggests there should be no disk either, which means the field 
is unlikely to be tangled on small scales. Further study of this source to search for evidence 
of outflow or a protostellar condensation could provide some support for the suggestion that 
this is a starless core. However, this would not completely explain the pattern, since all the 
starless cores observed thus far with SCUBA do exhibit depolarization toward their intensity 
peaks (Ward-Thompson et al. 2000). The cores observed by Ward-Thompson et al. (2000) 
are all significantly closer than Bl (140-170 pc). The total flux from one of these cores (L183) 
is 2.8 Jy at 800 /im (Ward-Thompson et al. 1994). Therefore, cores of comparable bright- 
ness, both closer and further (i.e., LBS 23's cores, Matthews, Fiege & Moriarty-Schieven 
(2002), hereafter Paper III) than Bl, do exhibit depolarization to the highest intensities. 
The p versus I relation of Bl-c is thus not easily explained purely by differing resolutions. 



4-1.2. Polarization in Individual Cores 

It is worth noting that the degree of alignment across the Bl-a, Bl-c, and Bl-d cores is 
particularly strong. Figure 3 shows that these three cores are not significantly asymmetric, 
so the polarization patterns cannot be said to align with any preferred axes of the cores. 
The elongated core Bl-b is the only one which exhibits systematic variation in polarization 
postion angle. Since this core is composed of two sources, the polarization patterns could be 
different within each core and then we observe their vector sum where the cores overlap. 

Several models exist which predict the polarization patterns across cores depending on 
various physical interpretations. A recent publication by Padoan et al. (2001) predicts the 
continuum polarization for protostellar cores assembled via supersonic magnetized turbulent 
flows in models of molecular clouds. They find that the universally observed trend of declin- 
ing polarization percentage with increasing intensity in star-forming cores can be reproduced 
by their model if grains are aligned only up to a threshold extinction (Ay ~ 3 mag in their 
simulation). Several of the resulting p versus / plots do suggest a flattening in the distribu- 
tion at high values of intensity. However, one difference between their simulations (see Figs. 
6 and 7 of Padoan et al. (2001)) and Figure 4 is the population of the bottom left quadrant 
(low p and low /); vectors exist in this region in the simulations, whereas this region is 
devoid of vectors in Figure 4, despite complete sampling down to the intensity threshold. 
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Measurement of low levels of polarized emission toward lower intensities will become possible 
with the next generation of detectors (e.g. SCUBA-2). 

The fact that there is a potential for confusion between two clouds along the line of 
sight makes interpretation of the polarization data more complicated in Perseus. However, 
constraints on the column density and spatial scale of potential variations in extinction of 
a second cloud suggest that the overwhelmingly predominant source of polarized emission 
is the Bl cloud and its associated cores. The models of Padoan et al. (2001) and Fiege 
& Pudritz (2000) relate to cores forming from lower density, filamentary structures. Bl is 
not a significantly elongated cloud (although the Perseus complex itself appears elongated 
on very large scales), and its cores appear to be forming from density enhancements in its 
cold interior. A recent model of magnetized cores predicts a relation between the geometry 
of a core and the measured polarization position angles (e.g. Basu (2000)). However, the Bl 
cores are not significantly elongated; the only asymmetry is in Bl-b, which is resolved into 
two sources with separate extended envelopes by Hirano et al. (1999). 



4.2. An Estimate of the Total Field Strength and Direction in Bl 

Myers & Goodman (1991) describe a method by which the total field strength and 
direction can be estimated in a cloud toward which a series of independent Zeeman measure- 
ments of Bi os and polarization measurements to infer the orientation of B pos have been made. 
Their formalism is described for absorption polarimetry, but is easily adapted to emission 
polarimetry if simple assumptions are made about the relation between the orientations of 
polarization vectors and the plane-of-sky magnetic field. 

A distribution of polarization data which can be fit by a single Gaussian is a good 
approximation to the precise function described in Myers & Goodman (1991) in the case 
where the nonuniform field is relatively small and the three dimensional random case can 
be assumed. As described in § 3.1 above, the polarization vectors associated with the faint 
lower density gas off the Bl cores can be fit with a single Gaussian distribution which has 
a mean of 91° (east of north) and a width of 19°. At 12" sampling, these vectors are not 
completely independent (since the JCMT beamwidth is 14"), but we use them to get a rough 
estimate of the field properties. 

The Myers & Goodman model is based on the presence of both a uniform (B ) and 
non-uniform (B n ) component to the field in the region. Non-uniform in this case refers 
to a disordered, possibly turbulent, field component (as opposed to an ordered, but not 
unidirectional field). Given a distribution of polarization data with mean polarization angle, 
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6, (in degrees) and dispersion, s, (in radians) and a series of Zeeman measurements with 
mean uniform component B 0z and dispersion a Bz1 the following quantities can be estimated: 
the inclination, 

% = arctan ( — ) ; (1) 



the total uniform mean field strength, 



sBn 



Bo = (2) 

cos 2 



and the dispersion in B, 

a B = a B N 1 ' 2 (3) 

where N is the number of correlation lengths of the non-uniform field component through 
the cloud. The correlation length is an expression of how quickly the non-uniform component 
of the field changes through the depth of the cloud. Values of B n separated by less than a 
correlation length are likely to be correlated, while those more spatially separated than a 
correlation length are independent. Myers & Goodman (1991) derive the relation: 

The maximum number of correlation lengths can be estimated under the assumption that the 
magnetic and gravitational energy densities are equal. Under this condition, Ay jB ss 0.17 
mag yuG -1 , where B is the total field strength (Chandrasekhar & Fermi 1953; Myers & 
Goodman 1988). This estimate yields N max pzs 12. 

The r.m.s. field strength and the relative strengths of the non-uniform to uniform mag- 
netic field energy densities can then be calculated according to the relations: 

< B 2 >V 2 = (B* + to*) 1 ' 2 ; (5) 

and 



B 



M n _ 3a 2 

where we have assumed a three-dimensional non-uniform field component. 



(6) 



We use the detection of — 27±4 fiG toward the position of IRAS 03301+3057 (Goodman 
et al. 1989) as an estimate of the B 0z value. In the absence of other independent detections, 
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there is no estimate of the dispersion in the line-of-sight field strength. Thus, we adopt a 
parameterized approach to combining the Zeeman and dust polarization data to estimate the 
three-dimensional field, where the ratio of <jb z /B 0z takes on a range of values (0.2, 0.4, 0.6 
and 0.8). Table 4 then summarizes properties of the magnetic field in Bl based on equations 
(1), (2), (3), (5), and (6) for iV = 1 and N = 10. 

Finally, the direction of the uniform component of the magnetic field can be estimated: 

Bq = (N cos 6_b + E sin Q B ) sin i + cos iz. (7) 

For emission polarization data, we assume the mean magnetic field direction, Ob, is related 
to G by ±90°. Emission local to each dust grain should be related to the local field direction 
in this manner (Hildebrand 1988), but for this to be the case in the vector averaged sum of 
all polarizations through the cloud (which is what we measure at the telescope) is to assume 
that the magnetic field orientation does not substantially vary through the depth of the 
cloud. This has been shown to be a poor assumption in some regions, where the polarization 
data support a more complex field geometry (i.e., OMC-3 in Orion A, Paper II; NGC 2024 
in Orion B, Paper III; and NGC 2068, Paper IV). However, Bl is a dark cloud and not part 
of a giant molecular cloud complex like Orion; therefore, a uniform field structure (at least 
away from the cores) is not an unreasonable first-order assumption here. We note however 
that in utilizing this relation in Bl, we are also assuming that all the polarized emission is 
arising in the Bl core, as opposed to in a second cloud as discussed above. 

Using equation (7) and a Bz /B 0z = 0.2 (hence % = 31° from the line of sight), 

B = 31//G ±(0.52JV + 0.0LE) -0.865 . (8) 

We note that the direction of the field along the line of sight is toward the observer. 

The field components can also be expressed in terms of two dimensions, x and z, where 
x lies along the plane-of-sky mean field direction, as estimated from the mean polarization 
direction. In this case, 

B Qx = B 0z tam (9) 

and substitution of B 0z and i gives: 

B = B 0x x + B J (10) 
= ±19 fiG x - 27 fxG z. (11) 



Based on comparisons with theoretical predictions based on the assumptions of equality 
of magnetic fields and kinetic energy and equivalence of magnetic and gravitational energies, 
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Crutcher et al. (1993) found that the £-field component measured in Bl agreed well with 
predicted values if the inclination of the field to the line of sight are close to zero. Further- 
more, because Bl was the only cloud with a detected magnetic field out of 12 in their survey, 
there was a concern that Bl might have an atypically strong field. Based on statistical anal- 
ysis, Crutcher et al. (1993) concluded that this need not be the case if the magnetic field 
in Bl lies nearly along the line of sight. The B 0z value of —27 ± 4 /iG may be substantially 
higher than the magnetic field strength at locations away from the dense cores. Therefore, 
our analysis represents upper limits to the field strenghths in low column density gas. 

However, our detection of ordered polarization vectors from dust associated with the 
main core of molecular gas in Bl indicates that at least some of the magnetic field in the 
region lies in the plane of the sky. Unfortunately, there is no way to unambiguously determine 
the plane-of-sky field strength from polarization data alone, since the degree of polarization 
may depend on grain size, shape, composition and degree of alignment or spin, as well as 
field strength (Hildebrand et al. 2000). We note that the OH Zeeman measurements can 
effectively probe regions with densities as high as n(i?2) ~ 10 4 cm" 3 (Crutcher et al. 1994), 
quite comparable to the densities associated with the dust emission. Assuming a temperature 
of 20 K for the main core, the 850 /iin flux density (at a level of 0.2 Jy) implies a column 
density of N(H 2 ) ~ 2 x 10 22 cm -2 , which corresponds to a density of n(i^2) ~ 3 x 10 4 cm -3 
if the emitting dust extends over the main core diameter of 0.8 pc (Bachiller et al. 1990). 
In reality, the emitting region of dust may be more confined along the line of sight, which 
would imply even higher densities. This is certainly the case where flux densities are high. 
The column density toward Bl-c is N(H 2 ) ~ 3 x 10 23 cm -2 . Since Zeeman data have been 
obtained toward regions of high density, it is possible that the OH Zeeman data and the 
dust emission polarimetry do not arise in precisely the same spatial regions of the cloud. 
They could, therefore, be sampling different field geometries, or at least different total field 
strengths. However, we have applied the Myers & Goodman (1991) method over a size 
scale comparable to the Arecibo beam at the frequency of the OH Zeeman measurements 
of Goodman et al. (1989), making our calculation quite reasonable for the lower column 
density dust. 



5. Summary 

We have detected polarized emission at 850 /zm arising from the dense interior of the 
Barnard 1 dark cloud. Our observations are centered on the molecular "main core" observed 
by Bachiller et al. (1990), in which three ammonia peaks were identified. Sub millimeter 
emission is detected coincident with each of the ammonia peaks. In total four dust cores are 
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identified, one of which has been resolved into two sources (Hirano et al. 1999). Two of the 
dust condensations, Bl-c and Bl-d, have not been previously observed in continuum. The 
Bl-a core is likely the 850 /xm counterpart of IRAS 03301+3057. This source appears quite 
faint at long wavelengths. We note that the detection of two new dense dust condensations, 
plus the Bl-b binary sources identified as young Class sources by Hirano et al. (1999), 
increases the number of YSOs or pre-protostellar objects in Bl by almost a factor of 2. This 
indicates that the star formation efficiency in Bl is likely much larger than the 0.5% estimate 
by Bachiller et al. (1990) since there may be other as yet unobserved pre-protostellar or 
protostellar cores in the cloud of which we have observed only a fraction. 

The polarized emission can be separated somewhat arbitrarily into two sub-groups by 
the coincident flux levels. Strong polarizations are measured toward faint dust emission 
regions where the mean polarization percentage is 4.5% (standard deviation 2.3%). The 
position angles are distributed about 90° (east of north) and can be fit by a Gaussian of 
mean 91.3° and dispersion 19°. The polarizations associated with high intensities (i.e., the 
cores) show smaller polarization percentages, with a mean of 2.6% and standard deviation 
1.4%. The vectors show alignment across the cores, but each core does not exhibit the same 
mean position angle. A comparison of each vector to adjacent values shows that vectors 
are strongly aligned with their neighbours in position angle. The largest discrepencies are 
observed at the "boundaries" between the dense cores and the lower column density dust 
emission in which they are embedded. 

Over the whole mapped area, we see evidence of the depolarization observed toward 
many star-forming cores. Interestingly, when the polarization percentages are plotted against 
intensity for individual cores, the p versus I relation flattens out at 30% and 40% of the peak 
emission from the Bl-c and Bl-b cores, respectively. The Bl-c core exhibits only one vector 
with a value of p < 1%, and thus the flattening in that core is real and not just an artifact 
of a lower limit on detectable polarizations or increasing optical depth. The observation 
of depolarization at the highest intensities of cores closer (Ward-Thompson et al. 2000) 
and further (Paper III) than Bl makes it unlikely that this effect is directly related to our 
resolution of Bl at 330 pc. 

None of the orientations of polarization vectors measured by SCUBA are directly re- 
lated to the two mean magnetic field directions detected with optical polarimetry of the 
Perseus complex (Goodman et al. 1990a). In the case where the Bl dense cores could be 
completely depolarized, the polarized emission along the line of sight to those cores would 
arise completely in the foreground cloud (proposed to be at 200 pc). If such a cloud contains 
fluctuations in extinction, those fluctuations must be on scales similar to the separation of 
cores in Bl at approximately half the distance. This is required to account for the differing 
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orientations measured in the two bright cores, Bl-b and Bl-c. We dismiss this scenario 
as unlikely since, unless the polarized emission from the foreground cloud rises in such a 
way as to offset the increasing intensity toward the Bl-c core peak, we should see a steeper 
depolarization toward Bl-c than in typical cores, not the threshold we observe. 

Finally, following the method of Myers & Goodman (1991), we have estimated the net 
field geometry in the Bl main core using our polarized emission data toward faint regions 
(centered on 90° east of north) and the line-of-sight field strength toward Bl measured by 
Goodman et al. (1989). We find that the total uniform field component is described by: 



B r 



31//G ±(0.52A> - OME) - 0.86i] , (12) 



under the assumption of a Bz /B 0z = 0.2. The ratio of the magnetic energy of the non- 
uniform component of the field to the uniform component ranges from 0.09 to 0.9 for this 
case, depending on the number of correlation lengths of the non-uniform component through 
the cloud. This result is roughly consistent with the theoretical predictions based on virial 
and magnetic equilibrium in the cloud, for which the line-of-sight field was comparable to 
the total predicted field. The high degree of ordering in the polarization data itself suggests 
that some component (possibly a significant amount) of the magnetic field could lie in the 
plane of the sky. 
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Table I. Observing Parameters for Jiggle Mapping of Bl 



Pointing 


; Center 


Number of 


R.A. (J2000) 


Dec. (J2000) 


Times Observed 


03 h 33 m 17!9 


+31°09 , 32"3 


16 


03 h 33 m 13:6 


+31°08'28'. / 25 


12 


03 h 33 m 18!3 


+31°07'03'. / 8 


30 



Note. - - The chop throw used for all observa- 
tions was 120" at a chop position angle of 65° (east 
of north). 
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Table 2. Barnard 1 850 /im Polarization Data 



A R.A. a 


A Dec. 


V 


dp 


Op 


e 


d6 


n 


(") 


(%) 


(%) 




n 


(°) 


Vectors associated with I > 720 mJy beam lb 




-58.5 


-136.5 


3.74 


0.56 


6.7 


-85.3 


4.3 


-64.5 


-136.5 


5.24 


0.83 


6.3 


-78.3 


4.5 


-52.5 


-130.5 


1.86 


0.50 


3.7 


68.1 


7.7 


-58.5 


-130.5 


2.36 


0.46 


5.2 


85.0 


5.6 


-64.5 


-130.5 


3.85 


0.59 


6.5 


-80.9 


4.4 


-70.5 


-130.5 


9.06 


0.76 


11.9 


-80.6 


2.4 


-46.5 


-124.5 


2.49 


0.63 


3.9 


67.3 


7.3 


-52.5 


-124.5 


3.10 


0.46 


6.7 


64.0 


4.3 


-58.5 


-124.5 


1.79 


0.44 


4.1 


85.4 


7.1 


-64.5 


-124.5 


2.29 


0.61 


3.7 


-72.8 


7.6 


-46.5 


-118.5 


4.10 


0.52 


7.9 


79.3 


3.6 


-52.5 


-118.5 


3.68 


0.41 


9.0 


79.8 


3.2 


-58.5 


-118.5 


2.25 


0.43 


5.2 


-87.8 


5.5 


1.5 


-112.5 


2.50 


0.51 


4.9 


-89.0 


5.9 


-4.5 


-112.5 


2.47 


0.46 


5.4 


-88.7 


5.3 


-10.5 


-112.5 


4.79 


0.43 


11.0 


66.7 


2.6 


-16.5 


-112.5 


4.52 


0.48 


9.5 


89.1 


3.0 


-22.5 


-112.5 


1.74 


0.50 


3.5 


85.6 


8.2 


-46.5 


-112.5 


4.25 


0.55 


7.7 


84.6 


3.7 


-52.5 


-112.5 


2.25 


0.47 


4.8 


78.3 


6.0 


-58.5 


-112.5 


1.46 


0.48 


3.0 


-66.1 


9.5 


-64.5 


-112.5 


2.13 


0.54 


3.9 


-80.1 


7.3 


13.5 


-106.5 


4.00 


0.50 


8.0 


89.8 


3.6 


7.5 


-106.5 


2.59 


0.35 


7.5 


89.2 


3.8 


1.5 


-106.5 


1.05 


0.32 


3.3 


81.9 


8.7 


-4.5 


-106.5 


1.57 


0.32 


4.8 


-82.1 


5.9 


-10.5 


-106.5 


2.75 


0.38 


7.2 


-86.7 


4.0 


-16.5 


-106.5 


4.26 


0.44 


9.7 


-86.8 


3.0 


-22.5 


-106.5 


3.99 


0.43 


9.3 


-89.1 


3.1 
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Table 2 — Continued 



A R.A. a 


A Dec. 




dp 




e 


dO 


\ j 


(") 


(%) 






(°) 




13.5 


-100.5 


1.28 


0.32 


4.0 


-80.5 


7.2 


-10.5 


-100.5 


2.08 


0.33 


6.4 


83.2 


4.5 


-16.5 


-100.5 


2.86 


0.41 


7.0 


-78.2 


4.1 


7.5 


-94.5 


1.13 


0.16 


7.2 


-63.3 


4.0 


1.5 


-94.5 


1.01 


0.15 


6.6 


-49.0 


4.4 


-10.5 


-94.5 


1.28 


0.36 


3.5 


47.2 


8.1 


-16.5 


-94.5 


2.53 


0.45 


5.6 


-82.8 


5.1 


19.5 


-88.5 


2.38 


0.41 


5.8 


33.8 


4.9 


-4.5 


-88.5 


1.18 


0.24 


5.0 


46.4 


5.8 


-10.5 


-88.5 


3.40 


0.41 


8.4 


58.5 


3.4 


13.5 


-82.5 


1.09 


0.28 


3.9 


28.0 


7.4 


7.5 


-82.5 


1.15 


0.18 


6.5 


58.9 


4.4 


1.5 


-82.5 


1.35 


0.18 


7.4 


67.5 


3.9 


-4.5 


-82.5 


1.52 


0.26 


5.8 


68.0 


4.9 


-10.5 


-82.5 


4.02 


0.46 


8.8 


66.9 


3.3 


13.5 


-76.5 


1.53 


0.37 


4.1 


-50.2 


6.9 


7.5 


-76.5 


1.50 


0.22 


6.8 


73.7 


4.2 


1.5 


-76.5 


2.35 


0.20 


11.9 


66.6 


2.4 


-4.5 


-76.5 


1.91 


0.27 


7.0 


71.2 


4.1 


-10.5 


-76.5 


4.48 


0.51 


8.8 


65.5 


3.2 


13.5 


-70.5 


3.00 


0.64 


4.7 


-80.1 


6.1 


7.5 


-70.5 


4.50 


0.37 


12.1 


-87.4 


2.4 


1.5 


-70.5 


1.98 


0.28 


7.2 


76.8 


4.0 


-4.5 


-70.5 


1.35 


0.34 


3.9 


64.7 


7.3 


7.5 


-64.5 


5.30 


0.65 


8.2 


87.8 


3.5 


1.5 


-64.5 


4.13 


0.45 


9.3 


77.1 


3.1 


-4.5 


-64.5 


2.98 


0.51 


5.9 


71.7 


4.9 


-28.5 


13.5 


3.56 


0.42 


8.5 


61.2 


3.4 


-34.5 


13.5 


1.25 


0.37 


3.4 


72.4 


8.4 


-22.5 


19.5 


2.04 


0.47 


4.4 


-56.5 


6.5 
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Table 2 — Continued 



A R.A. a 


A Dec. 




dp 


(Jr.-, 


e 


d9 


H 


\ j 


(%) 






(°) 




-28.5 


19.5 


1.22 


0.33 


3.7 


9.7 


7.8 


-34.5 


19.5 


1.41 


0.25 


5.6 


27.9 


5.1 


-40.5 


19.5 


1.45 


0.22 


6.5 


19.8 


4.4 


-22.5 


25.5 


3.12 


0.43 


7.2 


-24.0 


4.0 


-28.5 


25.5 


2.62 


0.31 


8.3 


24.8 


3.4 


-34.5 


25.5 


1.77 


0.19 


9.1 


31.6 


3.1 


-40.5 


25.5 


1.69 


0.17 


9.9 


31.9 


2.9 


-28.5 


31.5 


1.83 


0.31 


5.8 


30.0 


4.9 


-34.5 


31.5 


1.34 


0.18 


7.4 


44.9 


3.9 


-40.5 


31.5 


1.86 


0.20 


9.4 


47.0 


3.0 


-46.5 


31.5 


1.88 


0.30 


6.3 


50.3 


4.5 


-52.5 


31.5 


4.29 


0.48 


9.0 


64.6 


3.2 


-34.5 


37.5 


2.15 


0.29 


7.4 


47.1 


3.9 


-40.5 


37.5 


1.77 


0.30 


5.9 


46.9 


4.9 


-46.5 


37.5 


1.58 


0.38 


4.1 


40.6 


6.9 


Vectors 


associated 


with I 


< 720 


mJy 


beam~ lb 




-10.5 


-139.5 


3.44 


0.81 


4.2 


47.8 


6.8 


-22.5 


-139.5 


4.58 


0.87 


5.3 


-78.7 


5.5 


-34.5 


-139.5 


11.66 


0.93 


12.6 


88.1 


2.3 


-46.5 


-139.5 


6.88 


0.52 


13.2 


-70.8 


2.2 


-58.5 


-139.5 


4.05 


0.38 


10.6 


-65.0 


2.7 


-70.5 


-139.5 


8.00 


0.64 


12.4 


-73.1 


2.3 


13.5 


-127.5 


8.54 


0.96 


8.9 


89.9 


3.2 


1.5 


-127.5 


9.42 


0.57 


16.5 


-81.3 


1.7 


-10.5 


-127.5 


3.99 


0.41 


9.7 


87.1 


3.0 


-22.5 


-127.5 


3.64 


0.44 


8.2 


73.3 


3.5 


-34.5 


-127.5 


8.27 


0.45 


18.5 


-86.2 


1.5 


-46.5 


-127.5 


2.94 


0.33 


8.9 


67.2 


3.2 


-82.5 


-127.5 


5.19 


0.68 


7.7 


-88.7 


3.7 


25.5 


-115.5 


5.15 


0.64 


8.0 


-87.1 


3.6 
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Table 2 — Continued 



A R.A. a 


A Dec. 




dp 




e 


dO 


\ j 


\ j 


(%) 






(°) 




13.5 


-115.5 


8.13 


0.43 


18.7 


82.0 


1.5 


1.5 


-115.5 


2.91 


0.30 


9.6 


-85.2 


3.0 


-22.5 


-115.5 


2.37 


0.28 


8.3 


84.0 


3.4 


-34.5 


-115.5 


1.18 


0.32 


3.7 


56.5 


7.7 


-82.5 


-115.5 


1.77 


0.58 


3.1 


86.0 


9.4 


25.5 


-103.5 


6.52 


0.46 


14.1 


74.8 


2.0 


-34.5 


-103.5 


4.63 


0.31 


15.1 


81.2 


1.9 


-46.5 


-103.5 


2.64 


0.30 


8.7 


-85.3 


3.3 


-70.5 


-103.5 


1.46 


0.42 


3.5 


-89.6 


8.2 


-82.5 


-103.5 


2.14 


0.63 


3.4 


87.1 


8.4 


-22.5 


-91.5 


3.11 


0.31 


10.2 


89.5 


2.8 


-34.5 


-91.5 


3.41 


0.36 


9.4 


-87.3 


3.0 


-46.5 


-91.5 


1.16 


0.34 


3.4 


84.4 


8.3 


-58.5 


-91.5 


2.03 


0.32 


6.2 


79.1 


4.6 


-70.5 


-91.5 


1.65 


0.46 


3.6 


22.2 


7.9 


-82.5 


-91.5 


9.54 


0.92 


10.4 


-17.0 


2.7 


25.5 


-79.5 


3.43 


0.62 


5.5 


88.8 


5.2 


-22.5 


-79.5 


3.92 


0.51 


7.7 


75.4 


3.7 


-34.5 


-79.5 


3.78 


0.46 


8.3 


-85.8 


3.4 


-46.5 


-79.5 


4.09 


0.36 


11.5 


80.9 


2.5 


-58.5 


-79.5 


1.67 


0.32 


5.3 


35.8 


5.5 


-70.5 


-79.5 


3.06 


0.43 


7.2 


54.1 


4.0 


-10.5 


-67.5 


3.89 


0.31 


12.7 


83.2 


2.3 


-22.5 


-67.5 


6.48 


0.62 


10.5 


83.5 


2.7 


-46.5 


-67.5 


2.37 


0.41 


5.7 


82.5 


5.0 


-58.5 


-67.5 


4.19 


0.33 


12.8 


57.9 


2.2 


-70.5 


-67.5 


3.29 


0.49 


6.7 


15.0 


4.3 


13.5 


-55.5 


7.64 


0.86 


8.9 


0.1 


3.2 


1.5 


-55.5 


2.23 


0.43 


5.2 


-76.8 


5.5 


-10.5 


-55.5 


2.71 


0.41 


6.6 


70.5 


4.4 
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Table 2 — Continued 



A R.A. a 


A Dec. 




dp 




e 


dO 


\ j 


\ j 


(%) 






(°) 




-22.5 


-55.5 


8.82 


0.72 


12.2 


-85.7 


2.3 


-46.5 


-55.5 


3.04 


0.98 


3.1 


56.5 


9.2 


-58.5 


-55.5 


4.43 


0.81 


5.4 


-72.3 


5.3 


1.5 


-43.5 


2.54 


0.72 


3.5 


74.5 


8.1 


-10.5 


-43.5 


2.53 


0.59 


4.3 


86.9 


6.7 


-22.5 


-43.5 


7.08 


0.90 


7.8 


78.9 


3.7 


1.5 


-31.5 


6.44 


0.70 


9.2 


-70.8 


3.1 


-10.5 


-31.5 


5.43 


0.59 


9.1 


-67.9 


3.1 


-22.5 


-31.5 


9.13 


0.94 


9.7 


90.0 


3.0 


1.5 


-19.5 


8.57 


0.79 


10.9 


-84.0 


2.6 


-10.5 


-19.5 


4.18 


0.68 


6.1 


-60.5 


4.7 


-22.5 


-19.5 


4.86 


0.79 


6.1 


61.7 


4.7 


-34.5 


-19.5 


2.73 


0.67 


4.1 


-84.6 


7.0 


-46.5 


-19.5 


4.25 


0.76 


5.6 


-78.6 


5.1 


13.5 


-7.5 


3.89 


0.78 


5.0 


-68.9 


5.7 


1.5 


-7.5 


3.95 


0.76 


5.2 


-80.3 


5.5 


-10.5 


-7.5 


3.63 


0.58 


6.3 


-89.6 


4.6 


-22.5 


-7.5 


4.78 


0.54 


8.9 


-79.0 


3.2 


-34.5 


-7.5 


4.93 


0.45 


11.0 


-88.6 


2.6 


-46.5 


-7.5 


4.13 


0.62 


6.7 


-79.5 


4.3 


1.5 


4.5 


4.15 


0.74 


5.6 


66.9 


5.1 


-10.5 


4.5 


4.05 


0.55 


7.3 


-74.2 


3.9 


-22.5 


4.5 


1.25 


0.36 


3.5 


-51.8 


8.2 


-34.5 


4.5 


1.86 


0.30 


6.2 


54.0 


4.6 


-58.5 


4.5 


7.14 


0.87 


8.2 


42.9 


3.5 


13.5 


16.5 


8.55 


0.93 


9.2 


-58.2 


3.1 


1.5 


16.5 


3.23 


0.65 


5.0 


-63.8 


5.8 


-10.5 


16.5 


4.10 


0.45 


9.1 


-59.4 


3.2 


-58.5 


16.5 


4.97 


0.63 


7.9 


-86.8 


3.6 


1.5 


28.5 


7.91 


0.73 


10.8 


-84.6 


2.6 
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Table 2 — Continued 



A R.A. a 


A Dec. 


V 


dp 


Op 


6 


d6 


{ ) 


K ) 


l/o J 


l/o ) 




(°\ 


{ ) 


-10.5 


28.5 


4.22 


0.47 


9.0 


-77.2 


3.2 


-58.5 


28.5 


3.48 


0.43 


8.0 


68.6 


3.6 


-10.5 


40.5 


5.84 


0.76 


7.7 


-61.6 


3.7 


-22.5 


40.5 


2.24 


0.39 


5.8 


-3.1 


5.0 


-58.5 


40.5 


4.96 


0.54 


9.1 


72.1 


3.1 


-34.5 


52.5 


2.08 


0.69 


3.0 


-71.1 


9.5 


-58.5 


52.5 


2.88 


0.72 


4.0 


-47.1 


7.2 



a Positional offsets are given from the J2000.0 co- 
ordinates a = 03 h 33 m 20!9 and 5 = +31°09 / 03'/7 
(a = 03 h 30 m 15^0 and 5 = +30°59'00'. / in B1950.0). 
Vectors are binned to 12" sampling below the cho- 
sen threshold in total intensity and 6" sampling 
above. The total intensity at each vector position 
exceeds 20% of the faintest compact peak, Bl-d. 
This minimizes the chances of systematic effects 
from chopping to a reference position, as discussed 
in Appendix A of Paper II. 

b Using a calibration factor of 480 Jy beam -1 V -1 . 



Table 3. Peak flux densities at 850 /zm 



Core 


R.A. 


Dec. 


Speak 


Notes 




(J2000) 


(J2000) 


(Jy beam -1 ) 




Bl-a 


03 h 33 m 16!4 


+31 o 07'51" 


0.68 


IRAS 03301+3057 (Bl-IRS) 


Bl-bS 


03 h 33 m 2O 


+31°07'28" 


2.5 


Hirano et al. (1999) 


Bl-c 


03 h 33 m 17:7 


+31°09'31" 


3.0 


first continuum detection 


Bl-d 


03 h 33 m 16 ! :2 


+31°06'49" 


1.1 


first continuum detection 
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Table 4. Magnetic Properties of Barnard 1 



o B * i b B N a B <B 2 >- X I 2 M n /M. 



(pG) 


(degrees) 






(//G) 


(/iG) 




5.4 


31 


31 


1 


5.4 


32.4 


0.09 








10 


17.1 


42.9 


0.91 


10.8 


50 


42 


1 


10.8 


46.0 


0.20 








10 


34.2 


72.6 


2.00 


16.2 


61 


56 


1 


16.2 


62.6 


0.25 








10 


51.3 


105 


2.50 


21.6 


68 


71 


1 


21.6 


80.3 


0.28 








10 


68.4 


138 


2.80 



a For all calculations, the following values were used: \B 
27 /zG; s = 0.33 radians. 

inclination is measured from the line of sight. 



31 



Right ascension (B1950) 
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28 24 3:33:20 16 12 



Right ascension (J2000) 

Fig. 1. — 850 jum polarization pattern toward the "main core" of Barnard 1 is overlaid on 
the Stokes' / map. The greyscale range is —la to 3a. Polarization data were sampled at 
3" and have been binned to 6" (approximately half the JCMT beamwidth of 14") on the 
bright cores where signal-to-noise (a p ) is high and 12" in the fainter regions. All vectors 
are associated with Stokes' / values greater than 20% the Bl-d peak flux, a p > 3, and an 
uncertainty in polarization percentage, dp, < 1%. Red vectors have polarization percentage, 
p, < 1% and were not included in any analysis of the polarization pattern. The vectors 
are accurate in position angle to better than 10°. The peaks Bl-a, Bl-bN and Bl-bS have 
been labeled according to the designations of Hirano et al. (1999). We have designated the 
other two peaks Bl-c and Bl-d. The position of IRAS 03301+3057 is marked in blue where 
the lines denote the extent of the uncertainty in the IRAS position. The mean polarization 
percentage of the plotted vectors is 3.6% with a standard deviation of 2.2%. Coordinates in 
J2000 and B1950 are shown. 
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Fig. 2. — Distributions of position angle in Bl are shown using histograms of the 6 measured 
toward faint (dashed) and bright (solid) regions separately. The threshold used to discrimi- 
nate between faint and bright regions was a Stokes' / value of 720 mJy beam -1 , as discussed 
in §3.2. The vectors associated with lower column densities exhibit a Gaussian distribution 
with a mean of 91° and dispersion of 19°. The distribution of vectors toward the cores is not 
Gaussian, with each core dominating a different part of the total histogram. The position 
angles are labeled with the core dominating each peak. The dotted lines show the emission 
polarization position angles which correspond to the peaks observed in the Goodman et 
al. (1990a) distribution of optical absorption polarimetry (i.e., they are rotated by 90°). 
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Right ascension 

Fig. 3. — Core positions and the variation between neighbouring vectors. The positions of the 
cores are shown with contours of the Stokes' / emission. Contours plotted are: 0.2 to 1.0 Jy 
beam -1 in steps of 0.2 Jy beam -1 , and 1.4 to 2.6 Jy beam -1 in steps of 0.4 Jy beam -1 . (These 
correspond to visual extinction magnitude contours of 20 to 100 in steps of 20 magnitudes, 
and 140 to 260 in steps of 40 magnitudes.) We have compared each polarization vector shown 
in Figure 1 to its nearest neighbours in eight directions to a maximum radial separation of 
17". The data were then smoothed to 12" sampling. The length of the vectors plotted is 
the mean change Ap in the smoothed grid point, while the vector orientations indicate the 
mean change in direction A8. Zero change in angle is indicated by 0° (east of north). These 
data reveal that within cores, and between them, the position angles are strongly consistent, 
while the large vectors offsets from 0° at the boundaries show where the changes in vector 
character occur. Coordinates are J2000. 
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Fig. 4. — Depolarization in regions of low and high column density. Polarization is plotted 
against intensity on a log-log scale for our two populations of vectors. Low column density 
data are plotted with crosses while circles illustrate vector magnitudes associated with the 
cores. The slopes derived from fits to the p versus / profiles yield the same power law index, 
7 (where p = AP), plotted as the slope on the logarithmic scale. We note, however, that 
despite the scatter in the plot, the vectors at higher intensities seem to flatten in p. This 
can be seen more clearly in Figure 5. 
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Fig. 5. — Depolarization toward individual cores. To better examine the change in polariza- 
tion percentage at high intensities, we plot log p versus log I for individual cores. Essentially, 
we take all vectors around each core above the intensity threshold discussed in Figure 1. By 
plotting the three bright cores separately, it is clear that Bl-b and Bl-c show higher values 
of p than expected given the slope of the relation below 1 Jy beam -1 . The threshold corre- 
sponds to ~ 30% of the peak of Bl-c and ~ 40% of Bl-b. In the text, we discuss possible 
systematic effects which could explain this flattening (i.e., the truncation of the data set 
below 1% and optical depth). We conclude that the threshold is a real effect for the Bl-c 
core. 



